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ABSTRACT 

We show that feedback from active galactic nuclei (AGN) plays an essential role in 
reproducing the down-sizing phenomena, namely: the colour-magnitude relation; spe¬ 
cific star formation rates; and the a enhancement of early type galaxies. In our AGN 
model, black holes originate from Population ill stars, in contrast to the merging sce¬ 
nario of previous works. In this paper, we show how the properties of present-day 
galaxies in cosmological chemo-hydrodynamical simulations change when we include 
our model for AGN feedback. Massive galaxies become redder, older, less massive, 
less compact, and show greater a enhancement than their counterparts without AGN. 
Since we reproduce the black hole mass and galaxy mass relation, smaller galaxies do 
not host a supermassive black hole and their star formation history is affected very 
little, but they can get external enrichment from nearby AGN depending on their en¬ 
vironment. Nonetheless, the metallicity change is negligible, and the mass-metallicity 
relations, which are mainly generated by supernova feedback at the first star burst, 
are preserved. 

Key words: black hole physics - galaxies: evolution - galaxies: formation - methods: 
numerical - galaxies: abundances 


1 INTRODUCTION 


The A Cold Dark Matter (ACDM) cosmology is the cur¬ 
rently accepted best explanation for the structure we see in 
the Universe, and its evolution across cosmic time. On large 
scales, structure formation is driven by the gravitational in¬ 
fluence of dark matter, leading to Mpc-scale filaments, walls, 
and voids by the present day. Galaxies form within these, 
and their small-scale evolution is driven primarily by bary- 
onic processes. In the ACDM paradigm, all structure forms 
hierarchically; small halos collapse first, and merge over time 
to form ever larger ones. This has been challenged obser- 
vationally, howeve r , by the down-sizin g phenomenon (e.g. , 
Cowie et al] 19961 : (juneau et al.l 120051 : iBundv et al.l 1200a : 


Stringer et al.l 20091 1, whereby the mass of the most rapidly 
star-forming galaxies decreases with time. Similar phenom¬ 
ena have als o been observed as the evolution in mass-to- 
light rati os dT reu et all 2005 1. and the enhancement of a 


elements dThomas et al.l 2005 ) in early-type galaxies. In hy- 


drodynamical simulations, many galaxies that form at high 
redshift merge to form the most massive galaxies today, and 
so these galaxies tend to have ol der populations o f stars than 
pres ent-day low-mass gala x ies ( Kobayashi ct al. (200lj); see 
also iDe Lucia et al.l d2006l f ; iFontanot et al.l ( 20091 1 for semi- 


analytic models). With AGN feedback included, the star for¬ 
mation rate (SFR) at present can be r educed in line with 
observations (iTavlor fe Kobavashfeoiih , which will help re¬ 
produce the observations of the downsizing phenomenon. 


The potential importance of super massive black 
holes (BHs) in galaxy evolution has been made clear 
by the discovery that black hole mass correlates 
with numer ous properties of its host galaxy, including : 
bulge mass llMagorrian et al. 19981: Marconi fe Huntl [20031 : 


lHaring fe Rix[ 200j^_ Sani_etal. _[2011 ); stellar velocity dis¬ 


persion dFerrare se fe Merrittl l20 0ol: Icebhardt et ail I 2 OO 0 I : 


|TVemaine_et_aL |2002|; Kormendy_J^Jic 2013): luminosity 
i Kormendv fe Richstonelll995l : iMarconi fe Hunt||2003| ); and 
Serisc index / stellar c oncentration Jc raham et al.l 1200 ll : 
iGraham fe Drived boo 1 ?! ; [s avorgnan et al.ll2013fl . These cor¬ 
relations have been interpreted as im plyin g that galaxie s 
and BHs coevolve (but see also iJahnke fe Macciol 120110 . 
Such observed relations are reproduced by simulations 

dCroton et al.ll2006l: iDubois et al .11201 (j: ITavlor fe Kobavashil 

l20lA . as well as a better reproduction of the cosmic 
SFR (e.g..lBooth & Schavell2009l: ITavlor & Kobavashill2014l: 

IVogelshergeretalJb^ir 


We proposed a n ew model for AGN feedb ack in cos¬ 
mological simulations (ITavlor fe Kobavashill20l3) . With our 
AGN model, BHs form earlier in our simulations (z ~ 16) 
and with a lower mass (~ 10 3 M@) than previous works 
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(~ 10 5 Mq). This means that the effects of AGN feedback 
are different, especially at high redshift. However, the peak 
in cosmic SFR density still occurs at z = 2, consistent with 
observations and other theoretical works. Indeed, the ma¬ 
jority of stars in galaxies form before feedback from AGN 
becomes dominant, and so the scaling relations of galaxies 
are only weakly affected by AGN feedback. We will illustrate 
this with the mass-metallicity relations (MZR), which is the 
most stringent constraint on galaxy evolution. The galaxy 
size-mass relation is more affected, with AGN feedback in¬ 
creasing galaxy effective ra dius, but is still consistent with 
observations (but see also lDubois et all (120131 ) who find that 
AGN feedback is necessary to match observations). 

In this paper, we examine the effects of our AGN feed¬ 
back model on the chemodynamical evolution of galaxies in 
a cosmological context. We compare the z = 0 properties 
of galaxies in two cosmological simulations, one with, and 
one without, AGN with observations, and look in particu¬ 
lar at down-sizing and galaxy scaling relations. In Section [2] 
we introduce the code used to perform our simulations, and 
briefly describe our prescriptions for stellar and AGN feed¬ 
back. We outline the simulations performed in Section [3] as 
well as describing in more detail three example galaxies that 
will be used to highlight the differences in the simulations. 
Our results are presented and discussed in Section [4] and we 
draw our conclusions in Section [5] 


2 MODEL 


We use the s moothed parti c le hy d rodynamics (SP H) code 
GADGET-3 (ISpringel et al.l 1200 ll ; ISpringell 120051 1 . which 
is fully adaptive with individual smoothing lengths and 
times teps, and uses an entropy- conserving formulation of 
SPH ((Springel fc Hernauistll2002l l. This is augmented with 
various physical processes relevant to the formation and 
evolution of galaxies: metal-d ependent radiative co oling, 
star formation, stellar feedback (IKobavashi et al1l2007l l from 
AGB stars, Type la, and collapse supernovae a nd hyper¬ 
novae, and BH physics dTavlor fe Kobavashil 20 141) , T he BH 
physics is described in detail in Taylor fe Kobavashil (l2014h ; 
we summarize the key points below. 


2.1 Black Hole Physics 

i) BHs are seeded with a mass M Bee a from gas par¬ 
ticles of primordial composition (i.e. Z — 0) that are 
denser than some defined critical density (i.e. p s > p c ). In 
iTavlor fc Kobavashil (l2014l l. we found that values of M seet j = 
10 2 ~ 3 h _1 Mg and p c = 0.1 h 2 mn cm~ 3 were needed in or¬ 
der to produce results in agreement with observational con¬ 
straints (specifically, cosmic SFR density, the Mbh-u rela¬ 
tion, and the galaxy size-mass relation). Such a seed mass 
is consistent with the idea that the supermassive black holes 
observed today are the progeny of Population ill star rem¬ 
nants, in keeping with the lack of observational evidence of 
ch emic al enr ichm ent f rom pair-instability supernovae (e.g., 
iKobavashi et al1l201lh . 

Due to the limited mass resolution, M see d is much 
smaller than the mass of gas particles, M gas . For this rea¬ 
son, and in keeping with prev ious works (e.g.. lSpringel et al.l 
120051 : [Booth feSchavel 120091 ). we store the black hole mass 


separately, and use M gas as the dynamical mass until Mbh > 
M gas . Additionally, due to the spatial resolution, it is possi¬ 
ble for BHs to migrate out of galaxies in a few timesteps due 
to numerical effects. To avoid this, we calculate the centre 
of mass of all particles within the smoothing length of a BH 
in each timestep, and reposition the BH to that location. 

ii) BHs grow via both gas accretion and mergers. Given 
the spatial resolution of our simulations, it is necessary to 
adopt a sub-resolution model for both the growth of, and 
feedback from, BHs. To estimate the gas accretion ra te, we 
use the Bondi-Hoyle formula (Bondi & Hovlcl 1194411 mul¬ 
tiplied by a constant factor a that accounts for the fi¬ 
nite r esolution of the simulations. In ITavlor fe Kobavashil 
(1201411 . we treated a as a free parameter, and chose a value 
of 1 from observational constraints. M acc co uld be esti ¬ 
mate d by more sop h isticated methods (e.g., I Power et al.l 
l201ll : lAngles- Alcazar et al.ll2013l) . but higher resolution sim¬ 
ulations are required than in this work. We assume that ac¬ 
cretion is Eddington-limited at all times, and that BHs have 
a radiative efficiency of 0.1 dShakura fe SunvaevHl97~3h . 

Black holes also grow via mergers with one another. A 
pair of BHs merge if their relative speed is less than the 
local sound speed, and their separation is smaller than the 
gravitational softening length of the black holes. We use the 
sound speed as a measure of the local velocity scale, which 
precludes black holes that are undergoing a quick flyb y from 
merging (ISpringel et al.ll2005l : lDi Matteo et al]l200a ). 

iii) The feedback energy that a BH produces is pro¬ 
portional to the accretion rate, of which a constant frac¬ 
tion ef couples to neighbouring gas particles. With our AGN 
model, we found that ef = 0.25 is req uired to reproduce ob¬ 
servations dTavlor fe Kobavashill2C)f41 , while previous works 

have used 0.05 < _ ef < 0.15 (e.g.. ISpringel et al.l 120051 : 

iDi Matteo et al.l 120081 : iBooth fe Schava 200E ). The energy 
is distributed isotropically, kernel-weighted, and in a purely 
thermal form to a con stant number o f n eigh bouring gas par¬ 
ticles, Nfb = 72 (see [Tavlor fe Kobavashil d2014l) for more 
discussion). 


2.2 Stellar Physics 

Photo-heating is given by a uniform and time- evolving UV 
background radiation dHaardt fe Madaul Il996l l. and radia¬ 
tive cooling is computed using a metallicity-dependent cool¬ 
ing function ((Sutherland fe Dopital Il993l). The star forma¬ 
tion criteria used are the same as in Katz] dl992l l: (1) con¬ 
verging flow, (V • v)i < 0; (2) rapid cooling, t coo i < td yn ; 
and (3) Jeans unstable gas, tdyn < Gound- The star forma¬ 
tion timescale is taken to be proportional to the dynami¬ 
cal timescale (t s f = —tdyn), where c* is a star formation 
timescale paramete r, which we set to 0.02 . We spawn new 
star particles as in IKobavashi et ~ah ( 2007 1. and follow the 
evolution of the star particles at every timestep. The ener¬ 
gies are distributed to the same number of neighbouring gas 
particles as for AGN feedback, weighted by an SPH kernel. 
The energies of mass loss and supernovae are distributed 
in purely thermal form, although a fraction of it could, in 
principle, be distributed in ki netic form as a velocity per¬ 
turbation to the gas particles dNavarro fe Whitll 19931) . 

The c hemic al enri chmen t is computed with the 
scheme of IKobavashi (l2004f) . which does not include 
the instantaneous recycling approximation. Different from 
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IKobavashi et al.1 d2007h ’s simulatio ns, we adopt the initial 
mass function (IMF) of stars from Kroujsa (2008), a nd th e 
updated neucleosynthesis yields of IKobavashi et aD ( 2011 1 
for 1 — 50 Mq stars. Note that th e mass l oss from low-mass 
stars has been included since iKobavashil (|2004l l. The effects 
of hypernovae are included with a metal -dependent hyper¬ 
nova fraction as in IKobavashi fc Nakasatol (l201ll l. The pro¬ 
genitor model of Type la supernovae (SNe la) is very im¬ 
portant for predicting [a/Fe] of galaxies because SNe la pro¬ 
duce more iron than a elements (O, Mg, Si, S, and Ca). Our 
SN la model is based on the single degener ate scenario wit h 
the metallic ity effects of white dwarf winds llKobavashi et all 
1 19981 . l2000l l , and the lifetime distrib u tion f unctions are cal¬ 
culated as in IKobavashi fc Nomotol d2009ll . These provide 
an excellent agreement with the ob served elemental a bun¬ 
danc es in the solar neighbourhood (IKobavashi fe Nakasatol 
B)- As a result, half of iron is produced by SNe la, and 
the rest by core-collapse supernovae. 


3 THE SIMULATIONS 

In this section, we describe the set-up and large-scale evolu¬ 
tion of our simulations, before then describing the properties 
of three galaxies that will be used for illustrative purposes 
in Section Q] 



Figure 2. Cosmic star formation rate density as a function of 
redshift for the simulations with (solid black line) and without 
(dashed red line) AGN feedback. Sources of observational data 
(grey points) are listed in Section 13.11 


3.1 Cosmic Evolution 

We present simulations run both with and without 
the AGN model described in Section 12.11 With AGN 
feedback included, the fiducial model parameters of 
iTavlor fc Kobavashi! (120141 ) are used, namely a = 1, ef = 
0.25, M S eed = 10 3 /i _1 Mq, and p c = 0.1ft 2 mHcm -3 . The 
simulations in this paper are run at the same resolution 
with identical cosmological initia l conditions. We e mploy 
a WMAP-9 ACDM cosmology dHinshaw et al.l[2013 ) with 
h = 0.70, = 0.28, 11a = 0.72, fi b = 0.046, and a 8 = 0.82. 

The initial conditions are in a periodic, comoving, cubic vol¬ 
ume 25 ft -1 Mpc on a side, are chos en to have a central clus - 
ter at 2 = 0, and are different from IKobavashi et al.1 (|2007l ). 
which corresponds to a typical field of the universe. They 
contain 240 3 particles of each of dark matter and gas, with 
masses Mdm = 7.3 x 10' Mq and M gas = 1.4 x 10 7 Mq. By 
2 = 0 there are approximately 3 per cent more particles in 
total when AGN feedback is not included due to the forma¬ 
tion of star particles, and the influence of AGN feedback. 
We use gravitational softening lengths of cdm = 2.25ft -1 
kpc and e gaa = 1.125ft -1 kpc for dark matter and gas re¬ 
spectively. The finite volume of the simulation means that 
we do not necessarily obtain a representative sample of the 
most massive galaxies, while the low-mass end may be af¬ 
fected by the limited resolution of the simulations. When 
comparing to observations in the subsequent sections, our 
simulated galaxies are complete for 9 < log M/Mq < 11.4, 
-25.6 <M k < -21, and -22.1 < M B < 4 (B - V) 0 - 19, 
although we compare our simulations to each other across a 
wider range since their initial conditions are the same. 

The redshift evolution of the simulation with AGN is 
shown in Fig. [T| as well as the gas metallicity in the simula¬ 
tion without AGN. In both cases, a rich filamentary struc¬ 
ture exists (top row) at high redshift, to which star (second 


row) and black hole formation is mostly confined. Without 
AGN, the first star forms at z ~ 14 in the region that will 
subsequently collapse to form the largest cluster in the sim¬ 
ulation box. As early as 2 = 7, supernova feedback heats 
(fourth row) and chemically enriches the gas. This enrich¬ 
ment enhances star formation by allowing the gas to cool 
more rapidly. The cosmic SFR then peaks at 2 ~ 2 — 3. 
Supernova feedback causes galactic winds to be produced, 
which cause the intergalactic medium (IGM) to become en¬ 
riched by interstellar medium (ISM; sixth row). 

With AGN feedback included, BHs start forming from 
2 > 20, and prevent stars from forming until 2 ~ 11.5. When 
these first stars do form, however, they enrich the surround¬ 
ing gas and inhibit further BH formation. As in the case 
without AGN, the cosmic SFR peaks at 2 ~ 2, coincident 
with _a broader peak in_the BH accretion rate (Fig. 2 of 
ITavlor fe Kobavashil (120141) ). The BH number density peaks 
earlier, at 2 ~ 5.5, after which mergers reduce their num¬ 
ber. AGN and supernova feedback together produce stronger 
galactic winds than supernova feedback alone, and a larger 
region of the IGM is influenced (fifth row). The largest con¬ 
centration of stars and BHs is seen in the central cluster of 
galaxies, with a present total mass of ~ 10 14 ft -1 MQ, and 
which hosts the most massive black hole in the simulation 
with a mass of ~ 2 x 10 8 ft - 1 Mq. 


We show in Fig. [5]the cosmic SFR density for the simu¬ 
lations with (solid black line) and without (dashed red line) 
AGN feedback. Observational d ata a re plotted in grey, an d 


are taken from 

Bouwens et al. (201ll). Karim et al. (201ll). 

Cucciati et al.l ( 

20121). Oesch et al.. (20121), Burgarella et al. 

120131), Gunawardhana et al. (2013 

). and Sobral et al. 

(20131), adiusted to the Krouoa IMF 1 

Bernardi et al. 201(f); 

see Tavlor & Kobavashi (2014j) for a discussion on the dust 


extinction and IMF dependence. Both simulations show a 
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Figure 1. Evolution of cosmological simulation in 25h~ 1 Mpc box with (top 5 rows) and without (final row) AGN feedback included. 
We show projected gas density on the upper row, rest frame V band luminosity and stellar metallicity on the second and third, gas 
temperature on the fourth, and gas metallicity on the final two rows. 
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Figure 4. Star formation rate as a function of time since the Big 
Bang for the three example galaxies A, B, and C (Section |.‘1.2[ i in 
the simulations with (solid black lines), and without (dashed red 
lines) AGN feedback. 

peak in SFR of similar magnitude at 2 ~ 2 — 3, in agreement 
with observations. At higher and lower redshifts, however, 
the inclusion of AGN significantly reduces the SFR, and the 
simulated data lie much closer to the observations. At high 
redshift, although the BHs have not yet grown to be super- 
massive, gas has not been very much enriched, and so it is 
unable to cool efficiently. By low redshift, the supermassive 
BHs are able to accrete a large amount of gas, and provide 
enough energy to surrounding gas to drive powerful galactic 
winds. 

3.2 Individual Galaxies 

Galaxies are id entified using a Friend o f Fri en ds (FoF) 
algori thm dSpringel et al.l l200ll ) as in iTavlor fe Kobavashil 
((2014|), which produces a catalogue of groups of particles 
whose mass exceeds 32Mdm- For the analysis of properties 
of a galaxy, only those groups containing at least 50 star 
particles are used. In the simulation with (without) AGN, 
1380 (1009) groups are found by the Friends of Friends al¬ 
gorithm, of which 1125 (742) are classified as galaxies, cor¬ 
responding to 10 9 < M, /Mq < 10 12 . These differences are 
caused by AGN feedback at high redshift; although weak, 
early feedback from BHs suppresses star formation and low- 
mass galaxies survive. With AGN feedback included, 957 
galaxies contain at least 1 BH. 

Galaxy stellar masses and effective radii are esti¬ 


mated by fitting a core-Sersic profile llGraham et al.l 120031 : 
iTruiillo et al.l 120041) to the surface density. Approximately 
80 per cent of galaxies are well fit by such a function, and 
a further 10 per cent are well fit by a standard Sersic func¬ 
tion. In both cases we restrict the Sersic index to n = 4; this 
gives a much better fit to our simulated galaxies than with 
a Sersic function and varying n as in observations. Galax¬ 
ies that are not fit by either function tend to be near our 
resolution limit (~ 10 9 JW@) or merging systems, and are 
excluded from our analysis. Galaxy luminosities are esti- 
mated by applying the si mple stellar population models of 
iBruzual fc Chariot (i2003l ) to star particles with a Chabrier 
IMF, which gives very similar results to the Kroupa IMF 
used in our simulations. 

I 11 Section 0] we will highlight the effects of AGN feed¬ 
back by looking more closely at 3 example galaxies, matched 
between the runs with and without AGN. These galaxies, 
which we refer to as A, B, and C, have stellar masses ~ 10 12 , 
~ 10 11 , ~ 10 9 Mq respectively, evolve in different environ¬ 
ments, and are influenced differently by their BHs. We show 
in Fig. 0] images of the three galaxies at 2 = 0 from the 
simulation with AGN feedback, as well as their locations in 
the cosmological box of our simulation, colour-coded by in¬ 
tegrated V band luminosity. The global properties of these 
galaxies are summarized in Table [I] 

Galaxy A is the most massive galaxy in the simula¬ 
tions, and sits at the centre of the cluster. It experiences 
a major merger at 2 ~ 2, and numerous minor mergers 
throughout its life. With AGN, it hosts the most massive 
black hole, and AGN-driven outflows are produced inter¬ 
mittently. AGN feedback heavily suppresses star formation 
in the centre of A, as shown in the top panel of Fig. [4j al¬ 
though a small amount of star formation is ongoing in the 
outskirts of the galaxy (> 20/i _1 kpc) at a = 0 due to mi¬ 
nor mergers. As a result, the mass- and luminosity-weighted 
stellar ages are older, and the stellar mass of A is reduced 
by a factor of about 2.5. The removal of gas by AGN-driven 
winds reduces both M 200 and -R 200 , as well as the present 
gas fraction of this galaxy. This also changes the distribu¬ 
tion of heavy elements in this galaxy; stellar metallicity is 
slightly lower, though this difference is small enough that 
the mass-metallicity relations are not affected (Section 00]), 
and stars show a-enhancement. 

Galaxy B is a large field galaxy away from the main 
cluster (see Fig. 0 . It undergoes no major mergers, and 
very few minor ones. With AGN, an outflow is produced at 
2 ~ 2, and by the present, star formation is fully quenched 
over the whole galaxy by AGN feedback (Fig.0 since, unlike 
in galaxy A, there is no regular gas supply. There is an 
enhancement in both gas metallicity and temperature in a 
wider region of the IGM around B when AGN feedback is 
included. Its stellar mass is diminished by just over half by 
2 = 0, and its gas mass is reduced from 1.5 x 10 11 M© to 
10 10 Mq, causing its gas fraction to be reduced from 0.13 to 
0.03. Both stellar and gas-phase metallicities are lower with 
AGN feedback, and stars are a-enhanced. 

Galaxy C lies on a filament near the main cluster (see 
Fig. 0 and is the smallest of the 3 galaxies. It experiences 
a single merger at high redshift, but otherwise evolves pas¬ 
sively. AGN feedback has almost no effect on it; from the 
bottom panel of Fig. 0 the star formation rate within this 
galaxy is similar for the two runs across cosmic time, and the 
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Figure 3. Integrated V band luminosity image of the simulation with AGN feedback. The zoom-ins are 400h — *kpc on a side, and show 
the locations of galaxies A (left), B (bottom right), and C (top right), as well as their immediate environment in the simulation box 
(25fe _1 Mpc). 


Table 1.2 = 0 properties of galaxies A, B, and C both with and without AGN feedback. (1) Estimated from core-Sersic fit. (2) Gas 
fraction within R e , calculated as f g = M g /(M g + M* + Mbh)- (3) V band luminosity weighted within R e . (4) SFR weighted within R e . 
Note that there is no star-forming gas in A when AGN is included. (5) Galaxy colour in central I h~ 1 kpc. 


Galaxy 

M* 

/i^Mq 

(1) 

M200 
h- 1 M 0 

#200 
h~ 1 kpc 

/g 

(2) 

<Age> 

Gyr 

(3) 

[M/H], 

(3) 

[0/H] g 

(4) 

[O/Fe], 

(3) 

(B - V) 0 

(5) 

log (M B ; 

A (AGN) 

6.8ell 

9.5el3 

743.2 

0.001 

9.71 

0.23 

- 

0.18 

0.92 

8.4 

A (No AGN) 

1.7el2 

9.9el3 

752.6 

0.03 

2.19 

0.35 

0.39 

-0.10 

0.41 

- 

B (AGN) 

1.4ell 

5.1el2 

279.1 

0.01 

10.09 

0.10 

0.03 

0.20 

0.92 

7.9 

B (No AGN) 

3.7ell 

6.0el2 

295.8 

0.05 

2.90 

0.33 

0.41 

-0.10 

0.37 

- 

C (AGN) 

5.4e9 

l.Oell 

76.5 

0.17 

4.79 

0.15 

0.49 

0.05 

0.72 

5.3 

C (No AGN) 

4.8e9 

l.Oell 

77.4 

0.22 

3.08 

0.28 

0.35 

0.01 

0.54 

- 


differences are negligible. However, the chemical enrichment 
of its gas is enhanced with AGN feedback due to external en¬ 
richment from the strong, enriched winds generated in more 
massive galaxies. 


4 RESULTS 

In this section, we present the scaling relations of simulated 
galaxies in our cosmological simulations, one with the AGN 
feedback model of Section P included, and one without. 
Where applicable, we highlight the effects of AGN feedback 
by showing how the properties of galaxies A, B, and C (in¬ 
troduced in Section P change. 

4.1 Luminosity and Stellar Mass Functions 

We show in Fig. [5] the I\ band galaxy luminosity functions 
for our simulations. In the top panel, our calculated lumi¬ 
nosity functions, both with (solid black line) and without 


(dashed red line) AGN f eedbac k, are compared to the ob¬ 
servational data of lCole et al.l (l200ll . grey points) and the 
Schechter function fit to the data of iHuang et al.l (120031 . 
grey dot-dashed line). For the completeness limit of our 
simulation —25.6 < Mk < —20, discussed in Section ixn 
both simulations show fairly good agreement with observa¬ 
tions, though there are more low luminosity galaxies with 
AGN feedback (Section 13.21) . and fewer high luminosity 
galaxies, compared to the simulation without AGN. Look¬ 
ing more closely, we see that around the observed value of 
the Schechter function parameter M £ ~ —24.3, both sim¬ 
ulations produce very similar results, but are slightly lower 
than observed. This may suggest that supernova, and pos¬ 
sibly also AGN, feedback is slightly too strong around this 
mass in our simulations. The number of the most luminous 
galaxies (Mk < —25) is most affected by AGN feedback, and 
is decreased by AGN feedback, but is still slightly higher 
than in observations. This suggests that a larger simula¬ 
tion box, and possibly stronger AGN feedback at this mass, 
is required (see Section 0 for more discussion). Note that 
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Figure 5. Luminosity function for the simulations with (solid 
black line) and without (dashed red line) AGN feedback included. 
Our calculated l uminosity f unctio ns are compared with the obser¬ 
vational data of I Cole et al.l (l200lL grey p o ints) and the Schechter 
function fit to the data of I Huang et al ■ feool grey dot-dashed 
line). 



M b 

Figure 7. Colour-magnitude relation for our simulated galaxies 
with (black stars) and without (red diamonds) AGN feedback. 
Colours are measured within the central lh~ 1 kpc. Arrows show 
the motion of galaxies A, B, and C through the colour-magnitude 
plane when AGN feedba ck is inclu ded. Obser vational data (grey 
triangles) are taken from 1 Bender et al.l <119931) . 
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the agreement is not as good as in semi-analytic models 
(e.g., iBower et alJl200d : iGuo et al]|2013l : ISchave et al .11201-il) 
where the luminosity function is used to calibrate the param¬ 
eters of the subgrid physics. In hydrodynamical simulations, 
baryon physics is modeled on a particle basis with a smaller 
number of parameters. 

Fig. 1 shows the galaxy stellar mass functions for our 
simulat ions. Both simulati ons s how good agreement with the 
data of lLi fe Whitel (120091 ) and lBaldrv et all (120121) for the 
completeness limit of our simulations 9 < log M/M© < 11.4 
(Section 13.111 . The difference between the simulations with 
and without AGN is not as pronounced at the liigh-mass 
end as for the luminosity function, but there are only a small 
number of galaxies with log M/Mq > 11.5. The excess of 
low-mass galaxies when AGN feedback is included is due 
to the moderate suppression of SF early in the galaxies’ 
lives (see Section 13.21) . Despite this, our mass functions show 
a better fit to the observational data tha n obtained from 
other hydrodynamical simulations fe.g..lOp penheimer et all 
120101 : IPuchwei n fe Spr ingell 2013|; iKhandai et al.l 2014 . see 
also Fig. 5 of lSchave et alJ ll2014l) ). 


Figure 6. Stellar mass function for the simulations with (solid 
black line) and without (dashe d red line) AGN feedback. Obser¬ 
vational data are from iLi fe White! ||2009L grey dot -dashed line) 
who use the GAMA survey. and lBaldrv et all ( 20121 . grey points) 
who use SDSS. 


4.2 Colour — Magnitude Relation 

The colours of simulated galaxies are calculated by ap ply¬ 
ing s tellar population synthesis spectra llBruzual fe Charlotl 
1200311 to star particles as a function of age and metallic- 
ity. Although the dust content can, in principal, beesti- 
mated for simulated galaxies (e.g., iDe Lucia et al.l 120041 4. 
the reddening due to dust is not included in this paper. 
For this reason, we use optical bands where the bimodality 
llKauffmann et ahll2003l) is less visible than for UV bands. 
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Age [Gyr] 


Figure 8. Histograms of stellar age in the largest galaxy (A) in 
the simulation with (solid black) and without (dashed red) AGN 
feedback. 



°g(z/z 0 ) 


Figure 9. Histograms of stellar metallicity in the largest galaxy 
(A) in the simulation with (solid black) and without (dashed red) 
AGN feedback. 


Fig. [7] shows central B — V colour against B band mag¬ 
nitude for our simulations; black stars with, and red dia¬ 
monds without, AGN feedback. Colour is measured within 
the central l/i -1 kpc for galaxies, for co mparison with the 
observational data of iBender et al.l h 19931 ). Note that inte¬ 
grated colours across the whole galaxy tend to be bluer due 
to their colour gradient. In both cases, a distinct red se¬ 
quence ((B — V)o > 0.8), blue cloud ((B — V)o < 0.8), 
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Figure 10. Mean V band luminosity-weighted (top panel) and 
mass-weighted (bottom panel) ages of galaxies as a function of 
stellar mass. Error bars the standard error on the mean. Black 
stars denote the simulation with AGN feedback, and the red di¬ 
amonds without. 


and green valley ((B — F)o ~ 0.8) are seen. The colours 
of massive galaxies become redder wit h AGN fee dback , and 
consistent with the observations of lBender et al.l (ll993l . grey 
triangles). At Mb > —18, the simulated galaxies are on the 
extrapolated observed trend, but the number of these galax¬ 
ies seems to be larger than observed, which may be affected 
by the resol ution of our simulations. This result is very sim¬ 
ilar to iKobavashil (l2005i )’s simulations of elliptical galaxies, 
although AGN feedback is not included there. 

The movement of galaxies A, B, and C through the 
colour-magnitude plane is shown by the arrows. Galaxy A 
becomes much redder and less luminous, with AGN. Due to 
the significant reduction in star formation at low redshift 
(as shown in the top panel of Fig. a, and the dominance 
of very old stars (>10 Gyr) in the age distribution (Fig. [Sj) , 
this galaxy moves onto the red sequence when AGN is in- 
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eluded. Most of the young stars form in the outer regions of 
the galaxy as a result of mergers with gas-rich satellites. Fig. 
[9] shows the metallicity distribution of this galaxy, which re¬ 
mains mostly unchanged with the inclusion of AGN, except 
for the most metal-rich bin. The mass-weighted metallicity 
is not changed, but the luminosity-weighted metallicity is 
slightly decreased with AGN, and does not help in making 
colours re dder. If star formation is enhanced by AGN (e.g., 
Ism3l20iah . the metallicity, and thus colour, could be larger. 

Galaxy B also sees a significant reduction in star for¬ 
mation at late times when AGN feedback is included, to the 
extent that, by z = 0, star formation has been shut off al¬ 
most entirely (Fig. [4]). Unlike in A, there are no significant 
mergers to provide gas, and so star formation almost totally 
stops once the BH has heated or removed the gas from the 
galaxy. This allows B to move onto the red sequence. In 
galaxy C, the BH grows to only 2 x 10 5 /i -1 Mq, and is not 
able to significantly affect star formation or colour. 

Fig. llOl provides another representation of the age differ¬ 
ence. More massive galaxies tend to be older, and in partic¬ 
ular those with log M»/Mq >11 are older when AGN feed¬ 
back is included. This difference at log M*/Mq > 11 is more 
pronounced for luminosity-weighted ages (top pan el ). This 
AGN effect is in keeping with the findings of ICroton et ali 
d2006l) . but their galaxies were systematically older than our 
mass-weighted ages (bottom panel), which may be due to 
their cosmic SFR peaking at z ~ 4. 

4.3 Mass—Metallicity Relation 

We show in the previous section that the stellar metallicity 
distribution in the most massive galaxy was not much af¬ 
fected by the inclusion of AGN (Fig. O. We now show, in 
the top panel of Fig. 1111 the stellar mass-metallicity rela¬ 
tion (MZR) for galaxies in the simulations with and with¬ 
out AGN. In order to compare with observations of optical 
absorption lines, we show the V band luminosity-weighted 
stellar metallicity measured within the central 15/i _1 kpc of 
each galaxy. We note that, due to the metallicity gradient 
within galaxies, slightly higher values of [M/H]* are obtained 
with a smaller radius, but this effect is much smaller than 
the scatter in the simulated relation. There is no significant 
change in the MZR between the two simulations. This is 
because AGN feedback quenches star formation most effi¬ 
ciently at low redshift, after the peak of star formation and 
chemical enrichment. The simulated galaxies are in excellent 
agreement with the observations of Thomas et al.l d2005l ). 
with [M/H]* ~ 0.2 at M* ~ 10 n MQ. The scatter of the 
simul ations is also c ompa rable to that of observations. Note 
that lGallazzi et al.l (l2005l) ’s MZR of SDSS galaxies gives sys¬ 
tematically lower metallicities at all masses, and even lower 
metallicities below A/* ~ 10 10 Mq. Although these obser¬ 
vations are subject to large scatter, most of our simulated 
galaxies lie above this MZR. 

In the bottom panel of Fig. 1111 we show the stellar 
mass-iron abundance relation. A stellar mass-iron abun¬ 
dance relation is seen in our simulations because, although 
there is less SN la enrichment in massive galaxies (as we 
show in the [a/Fe] ratios of the following section), half of 
iron is produced from core-collapse supernovae, and thus 
the iron abundance also increases for massive galaxies. We 
see that this relation mirrors closely the MZR in the up- 
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Figure 11. Stellar mass metallicity relation for simulations with 
(black stars) and without (red diamonds) AGN feedback. Top 
panel: V band luminosity weighte d total metallicity in 1 5/2. 1 kpc. 
Observational data are taken from Th omas et ali j 20051 , solid blue 
and dashed green lines), and lGallazzi et al.l (12003. dot-dashed yel¬ 
low line; scatter indicated by dotted lines). Bottom panel: V band 
luminosity-wei ghted iron abundan ces in 1 kpc. The observed 
relation of ljohansson et al.l (12012! ) is also shown by the solid line. 
Arrows show the movement of galaxies A, B, and C through the 
MZR when AGN feedback is included. 


per panel, albeit with greater scatter. Howeve r, there is no 
[Fe/H],-M* relation in lJohansson et ahl j2012l Vs analysis of 
SDSS data, which may be due to the uncertainties in mea¬ 
suring elemental abundances from integrated spectra. 

The s lope of the MZR of our simulations is almost the 
same as in lKobavashi et ahl (120071) . but there is a systematic 
offset. This is because the Krou pa IMF is used i n this paper, 
rather than the Salpeter IMF in lKobavashi et ahl (120071) . Al¬ 
though the slope of massive stars is very similar between 
these IMFs, the difference for low-mass stars changes the 
norma lization , which ch ang es the net yields from 0.6Zq to 
1 Zq (IKobavashi et al]|201ll ). The slope of the MZR is gener- 
ated by the mass-dependent galactic winds (iKobavashi et al .1 
120071) . which is also the case in this paper. 

Arrows denote the movement of galaxies A, B, and C 
through the MZR. Galaxies A and B, in which star formation 
is efficiently quenched, move down the observed relation to 
lower stellar mass and slightly lower (~ O.ldex) metallicity. 
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Figure 12. Gas phase mass-metallicity relation. SFR-weighted 
gas metallicity for the simulations with (black stars) and with¬ 
out (red diamonds) AGN feedb ack. The solid line show s the 
SPSS observations, tak en from iKewlev fc Ellisonl <1 20081) with 
iKewlev & Dooi ta (2002)’s method, and arrows show how galaxies 
B and C move through the MZR when AGN feedback is included 
(there is no star-forming gas in A at z = 0 when AGN feedback 
is included). All metallicities are measured within 15h —1 kpc of 
the galactic centre. 


Figure 13. [a/Fe]—<7 relation for our simulations with (black 
stars) and without (red diamonds) AGN feedback. Arrows de¬ 
note the movement of galaxies A, B, and C through the [a/Fe]-cr 
relation when AGN fee dback is included . Observationa l data from 
I Kuntschner et al ] 1I2 OIOI, solid blue line ).IS nola.or et alJ d201(]l . dot¬ 
ted green line), and I Johansson et alJ ll2012l . orange dashed line) 
are also shown. V-band luminosity weighted [O/Fe] is measured 
within 1 effective radius for the simulated galaxies. 


Galaxy C, on which the effect of AGN feedback is minimal, 
moves only very slightly in either stellar mass or metallicity. 

We show also, in Fig. 1121 the gas-phase MZR, using 
SFR-weighted oxygen abundances for comparison with ob¬ 
servations that are estimated from emission lines. We again 
compare the simulations with and without AGN, and in¬ 
clude all of the gas in the central 15ft, _1 kpc of each galaxy. 
As with the stellar MZR, there is no significant difference be¬ 
tween these simulations. The simulated g alaxies are consis¬ 
tent w ith SP SS observat i ons, taken from IKewlev fc Ellisonl 
(I2OO8I) with IKewlev &; Dopital ll2002l) ’s method (shown as 
the solid line), but there is much greater scatter than ob¬ 
served. The scatter is als o greater than ha s been found 
in other simulations (e.g. IPave et al. llmH), though the 
adopted models differ both in their nucleosynthesis yields 
and the fact that our model includes feedback from AGN. 

The movement of B and C through the gas-phase MZR 
is within the scatter, and is more complicated. As well as in¬ 
fluencing the enrichment history of galaxies, AGN feedback 
also drives galactic winds from the most massive galaxies. 
This alters the metallicity of both the ISM and IGM com¬ 
pared to the case without AGN (Fig. [T]). Consequently, there 
is no trend on the movement of our example galaxies (note 
that there is no star-forming gas in galaxy A at z = 0, and 
so it has no arrow in Fig. m. 



Figure 14. Histograms of [a/Fe] in the largest galaxy (A) in the 
simulation with (solid black line) and without (dashed red line) 
AGN feedback. 
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4.4 [a/Fe] - Velocity Dispersion Relation 


The ratio of a-elements to iron abundances gives informa- 
tion abo ut the star formati on timescale of a galaxy (e.g., 
iKobavashi fe Nomotol 120091 1. Q-elements are mostly pro¬ 
duced by core-collapse supernovae with a short timescale 
(~Myr), while iron-peak elements are produced more by 
SNe la on longer timescales (~Gyr). Therefore, at the early 
stages of galaxy formation, only core-collapse supernovae 
contribute to the chemical enrichment of the galaxy, lead¬ 
ing to a higher a/Fe ratio relative to the solar abundance. 
At later stages, due to the delayed iron production by SNe 
la, [a/Fe] decreases as a function of time. The timescale of 
SNe la is an important issue, and depends on the progenitor 
model of SNe la; we take the best theoretical model that 
reproduces the observed [a/Fe]-[Fe/H] relation in the Milky 
Way Galaxy ('Section 12.211 . 

Fig. 1131 shows the [a/Fe]-<r relation for galaxies in our 
simulations with and without AGN feedback. [a/Fe] is mea¬ 
sured, V band luminosity-weighted, within the central effec¬ 
tive radius of a galaxy. Velocity dispersion, a, is also mea¬ 
sured within lR e of the centre. Without AGN, star forma¬ 
tion takes place longer in more massive galaxies because of 
their deep potential well, and thus [a/Fe] is lower for more 
massive galaxies, which is totally opposite to the observe d 
relations. This was also the case in IKobavashi et all (l2007tl . 
where AGN was not included. With AGN, it is possible to 
quench the recent star formation in massive galaxies. Al¬ 
though the effects of AGN feedback depend on environment, 
on average AGN feedback is more effective at quenching star 
formation in more massive galaxies. These galaxies are made 
up of older stellar populations than their counterparts with¬ 
out AGN, as can be seen in Fig. 1101 land is shown explicitly 
for the most massive galaxy in Fig. Since star forma¬ 
tion is most suppressed at low redshift, the majority of stars 
formed from gas that had not been as enriched with iron as 
in the case without AGN. As a result, more massive galaxies 
tend to have higher values of [a/Fe] when BHs are included; 
this is shown explicitly for galaxy A in Fig. 1141 The move¬ 
ment of galaxies A and B, in which AGN feedback is strong, 
in the [a/Fe]-cr plane also illustrate this, while C moves very 
little, since its star formation is not affected much. 


Observational data show a tight relation between [q/Fe] 
and g alaxy mass for b o th ce ntral values (iThomas et al.l 


I 2 OO 3 I: Jjohanssonet_ah| ]2012| ) and value s within R , 


( Kuntschner et al.ll201C ; ISpolaor et al.l l20ld ) , although the 
slope and absolute values depend on the analysis of the ob¬ 
served absorption lines. The aperture effect should be negli¬ 
gible as the observed radial gradient of [a/Fe] is very small. 

When AGN feedback is included, a number of 
galaxies move towards the observed relations for early- 
type galaxies dKuntschner et alJ l20ld : ISpolaor et alJ l20ld : 
Ijohansson et al] ~ 20121 . note that all simulated galaxies 
are displayed). For logo- > 2.1, the relation is reversed 
when AGN is included, more in keeping with observa¬ 
tions. However, the scatter of the simulations is still larger 
than observed. An improved model of AGN feedback may 
be re quire d , or o ther possibilities have been discussed in 
IKobavashi] (|2005f ). including (i) IMF, (ii) nucleosynthesis 
yields, (iii) binary fraction, and (iv) selective mass loss. 

The [a/Fe]-o- relation was als o shown with semi- 
analytic models with AGN feedback in lArrigoni et al] (l2010l l 
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Figure 15. Specific star formation rate against stellar mass for 
the simulations with (black stars) and without (red diamonds) 
AGN feedback. Arrows show the effect of AGN feedb ack on galax¬ 
ies A , B, and C. The fit to observational data by ISali m et al.l 
(1200711 is shown by the solid black line. 



andjGargiulo_et_aD (2014), both with a top-hea vy IMF, and 


Calura fe Mencil (2009) with a variable IMF. lYates et al.l 
(|2013h also requires modification of the SN la model. How¬ 
ever, our results suggest that the [a/Fe]-<r relation can ba¬ 
sically be reproduced with our AGN feedback and standard 
stellar physics (IMF, binary fraction, and SN la model). 


4.5 Specific Star Formation Rate 

Finally, in Fig. [15] we show the relationship between specific 
SFR (sSFR) and stellar mass for our simulated galaxies. In 
both simulations, we see a trend for more massive galax¬ 
ies to have lower sSFR, which is caused by supernova feed- 
back. This is in qualitative agreement w ith observations (e.g. 
ISalim et alfeOOTtlCiambur et al.l2013li . With AGN feedback 
included, the most massive galaxies, such as A and B, show 
significant reduction in sSFR. Smaller galaxies are less af¬ 
fected - for example, C moves very little in the sSFR-M* 
plane. 

Bo th obs e rvatio ns and semi-analytic modeling 
(ICiambur et al] l2013fl find a bimodal distribution of 
sSFR at z = 0, corresponding to low-mass galaxies with 
ongoing star formation, and quenched massive galaxies. 
Our results do not show such a distinction, which may be 
because there are only a small number of galaxies with 
log M*/Mq > 11 in our simulation box, and because the 
sSFR of < 10 -lo yr -1 in less-massive galaxies (~ 10 10 Mq) 
are affected by the resolution of our simulations (linear 
features in Fig. nsj. 
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5 CONCLUSIONS 

We have shown the effects of feedback from AGN on galaxy 
evolution using two cosmological hydrodynamical simula¬ 
tions, whose only difference is whether AGN feedback is in¬ 
cluded. In our AGN model, black holes originate from Pop¬ 
ulation ill stars, in contrast to the merging origin of previ¬ 
ous works, and grow via mergers and gas accretion. Cosmic 
SFRs are greatly reduced at both high redshift (z > 3) and 
low redshift (z < 1; Fig. 0 with AGN feedback. Since the 
black hole mass-galaxy mass relation is reproduced, mas¬ 
sive galaxies host massive BHs. These produce AGN-driven 
winds and quench the star formation there, causing them 
to be less massive, less compact, and redder by z = 0 than 
their counterparts without AGN. On the other hand, smaller 
galaxies do not host a supermassive BH and their star for¬ 
mation history is affected very little, but can get external 
enrichment from nearby AGN, depending on environment. 
Nonetheless, galaxy mass-metallicity relations are not af¬ 
fected much since most stars form before AGN feedback be¬ 
comes dominant. The [a/Fe] ratio of stellar populations is 
higher for massive galaxies due to their shorter star forma¬ 
tion timescale, without modifying stellar physics (i.e., IMF, 
nucleosynthesis yields, SN la progenitors, and binary frac¬ 
tion). Feedback from AGN plays an essential role in galaxy 
formation, namely in reproducing downsizing phenomena 
such as the colour-magnitude relation (Fig. 0, specific star 
formation rates (Fig. 1151) . and the a enhancement in early 
type galaxies (Fig. fTTTI) . 

Although our simulated galaxies are in better agreement 
with observations, there are still some mismatches seen, in¬ 
cluding: the luminosity function overpredicts the number of 
the most luminous galaxies; and the [a/Fe]-<r relation shows 
a larger scatter than observed. These suggest stronger feed¬ 
back to suppress star formation at low redshift (z < 0.5; 
Fig. 0. This cannot be achieved simply by changing the 
param eters o f BH phy sics i n our m odel, as these were cho¬ 
sen in iTavIor fe Kobavashil (l2014l i to provide the best re¬ 
production of observational constraints: the cosmic SFR; 
Mbh~<7 relation; and galaxy size-mass relation. There may 
be room for improvement by simultaneously changing the 
numbe r of feed bac k neig hbo ur pa rticles, as discussed in de¬ 
tail in lTavlor fe Kobavashil (2014). Alternatively, including 
smaller scale physics of AGN, kinetic feedback from jets, 
and induced star formation may be required, which will be 
investigated in future works. 
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